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The low temperature magnetization process of the ferromagnetic-antiferromagnetic Heisenberg 
chain is studied using the interacting boson approximation. In the low field regime and near 
the saturation field, the spin wave excitations are approximated by the 5 function boson gas for 
which the Bethe ansatz solution is available. The finite temperature properties are calculated by 
solving the integral equation numerically. The comparison is made with Monte Carlo calculation 
and the limit of the applicability of the present approximation is discussed. 
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§1. Introduction 

The spin-1/2 ferromagnetic-antiferromagnetic alter- 
nating Heisenberg chain has been int(^ai3i(i'v_s±ufiied re- 
lated with the Haldane gap problemMBBy yBLI'El' In the 
limit of strong ferromagnetic bond, this model tends to 
the spin-1 antiferromagnetic Heisenberg chain(AFHC). 
On the other hand, it reduces to an assembly of trivial 
isolated dimers if the ferromagnetic interaction vanishes. 
It also tends to the dimerized spin-1/2 AFHC which de- 
scribes the spin-Peierls state if the ferromagnetic bonds 
are replaced by the antiferromagnetic bonds. Thus the 
ground state of this model interpolates the Haldane gap 
state and the spin-Peierls state continuously. 

The magnetization process of this modeLLas ibpjsn also 
studied theoretically and experimentally .QU'Q'Q'B' It is 
pointed out that the critical behaviors at the lower criti- 
cal field and the saturation field are described by the free 
fermion theory, which, is equivalent to the boson with 
hard core interaction.^ It is the purpose of the present 
work to investigate the thermal effect on the magnetiza- 
tion process of this model. 

The Hamiltonian of the present model is given by 

N N 2N 

H = jY,S2iS2i+i + J'Y. '^2i-i52i - gmHY^Sf, 

1=1 1=1 1=1 

(1-1) 

where Si is the spin operator with magnitude 1/2 on the 
/-th site. It is assumed that J is positive (antiferromag- 
netic) and J' is negative (ferromagnetic). The external 
magnetic field, the electronic g factor and the Bohr mag- 
neton are denoted by g and , respectively. In the 
following, we take the unit J — 1 and g^-Q ^ 1- The 
ratio J'/J is denoted by /3. In the next section, the 
magnetization process of this model at low temperature 
T << 1 is studied by means of the mapping onto the 



interacting boson system which is solvable by the Bethe 
Ansatz method. The quantum Monte Carlo results are 
presented in §3. In the last section, the both results are 
compared and the limit of applicability of the interacting 
boson approximation is discussed. 

§2. Interacting Boson Approximation 

2. 1 Near the saturation field 

In the vicinity of the saturation field, we may apply 
the spin wave method starting from the fully polarized 
state. In order to incorporate the nonlinear terms in a 
convepisnt way, we apply the Dyson-Maleev transforma- 
tioiJaMEJ as follows, 

St = {\- a\ai), S+ = (1 - ajai)ai, S[- = aj, (2.1) 

to obtain the Boson representation of the Hamiltonian 
as 

H^Hq+H^j^^, (2.2) 

Uq = ^{- - -(4, - a2/+i)(a2/ - a2i+i)} 
1=1 

+ ^^E^i - 2^^21-1 - 4)(a2/-i - a2/)},(2.3) 
1=1 

1 ^ 

"^int = ^2 5I"2;«2/+i(«2i - a2i+if 
1=1 

- 2 X! 4i-i4/(a2;-i - a2if. (2.4) 

1=1 

By the Fourier transformation 

a2i = ^ Ck exp(ifcO, 
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-i = ;^E^fc^^P(i^(^ + ^))' (2.5) 



where the length of the unit cell (=twice lattice constant) 
is set equal to unity, Hq is transformed as, 



+ i(J(fc)4cfe + J*(fc)ct.4)|, 



dk) 



(2.6) 



where 

E,^il±p^-NH, J(fc)=ei^/2+/3e-i^/^ (2.7) 

After the Bogoliubov transformation which diagonalizes 
Ho, we have 



Ho = Eo + ^{S+(fc)4afe + E.{k)l3ll3k}, (2.8) 
fe 

where 

E±{k) ^-^(l + (3± v/l + /32 + 2/3cosfc) + i7, (2.9) 

and afc and f3k are boson operators given by the linear 
combinations of Ck and dk ■ The magnetization per unit 
cell M is expressed as 



M 



1 - ^ E{< "1"'= > + < >}■ (2.10) 



where <> represents the thermal average. The num- 
ber of the quasiparticles corresponds to the number of 
inverted spins. For H > 1, the excitation energy of 
the quasiparticles are positive corresponding to the fully 
magnetized ground state, while ior H < 1, the quasi- 
particles with momenta fc ~ tt have negative energy and 
these quasiparticles are present even at zero tempera- 
ture leading to the deviation from the saturated magne- 
tization. Therefore the saturation field Hg is given by 
Hg = 1. Near the saturation field at low temperatures, 
the density of the quasiparticles is low and they are dis- 
tributed around the minimum of the excitation spectrum 
at A: '--^ TT of -E+ . Therefore we only take into account the 
-|- branch and shift the origin of the momentum by n 
as k = n + q. In the following, we use g's instead of 
fc's as momentum indices. Expanding around q = and 
keeping the important terms at fc ~ tt, we find 



H = Eo+J2{-^ 



H 



-}alaq 



1 

2N 



91 ,92,93:94 



'^91 '^92 '^93 '^94 ^91 +92, 93 + 94- (2-11) 



Here the same Fourier and Bogoliubov transformations 
are applied to the interaction Hamiltonian Tijj^^ and the 
terms containing the — branch are neglected. The co- 
efficient of Hjjj^ is independent of /? at fc '--^ tt as far 
as < 0. This is due to the fact that the -I- branch 
represents the in-phase motion of the spins connected 
by the ferromagnetic bonds. The correlation length ^ 
is given by C = v^l /3 | /2(1+ | /3 |). The Hamiltonian 



(2.11) yields the one for the 5" = 1 AFHci) in the hmit 
/3 ^ — oo. This is the Hamiltonian of the (5- function bo- 
son which is [Solved by the Bethe Anzatz method, by Lieb 
and Liniger.O According to Yang and Yang,!!^ the fi- 
nite temperature properties of this model is described in 
terms of the solution of the integral equations. 



-A + q' 



tc 



& + 



^ln(l + e-^«)/*), 
(g- g') 



27rp(g)(l + exp(e(g)/t)) = 1 + 2c 



(2.12) 

P(g0rfg' 

c2 + (g -(?')"' 
(2.13) 



where 



1 



2T 



The magnetization per unit cell is given by 

/oo 
p{q)dq. (2.15) 
-oo 

This set of integral equations are solved numerically to 
obtain the magnetization curves. 

2.2 The low field regime 

In the low field regime < H ^ A = He, where 
A is the energy gap, the excitations from the ground 
state are assumed to be the spin-1 bosons with parabolic 
dispersion above the gap A at fc = TTj-and short range 
repulsion similarly to the spin-1 AFHCB^ Again shifting 
the momentum as k = tt + q, the effective Hamiltonian 
is given by 



2„2 



eq 



{hs- l)}aj „a 



- E 



t t 



9l,92,93,94,s,s' 
X Ctq3,s'Oiq4,s^qi+q2,q3+qiJ 



(2.16) 



where h = H/Hc- The suffices s and s' denote the z- 
component of the spin of the excitations and take the 
values 1, and —1. The values of A and ^ are estimated 
by extrapolating the exact diagonalizatiDn datEU'cP to 
iV — > oo by the Shanks' transformationliJ-' as shown in 
Fig. 1. Our boson Hamiltonian again yields the carre- 
sponding boson Hamiltonian for the S — \ AFHCh^ in 
the limit (3 — )■ — oo. 

The magnetization per unit cell is given by 



1 



9 



-1"9, 



-!>}. (2.17) 



We further neglect the interaction Cs^s' for s ^ s' . This 
is justified as follows. In the limit h « 1, and T « A, 
the occupation number of each bosonic state is low and 
all interactions may be neglected. On the other hand, in 
the vicinity of the lower critical field (ft, ~ 1), the occupa- 
tion numbers of the low energy states of s = 1 excitation 
are not necessarily small. Therefore the interaction ci^i 
must not be neglected. However, other interactions can 
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be neglected because the occupation numbers of other 
excitations are small at low temperatures. Nevertheless 
we keep the interaction c_i._i for the sake of symmetry 
between the excitations with s — 1 and s = — 1 . Follow- 
ing Takahashi and Sakaijj-' the value of the interaction 
constant Cs.s is identified as Cs.s = 1/2 by comparing the 
ground state magnetization for Hg » H » He ob- 
tained from the bosonic theory with the classically cal- 
culated magnetization M = H. Thus both s — 1 and 
s — —1 excitations are described as (5- function bosons 
and the magnetization curve at finite temperature is ob- 
tained using the nume rical solution of the integral equa- 
tions (Ell) and (^.131) with 



A = 



2{hs - 1) 



1 



IT 

<=— , (2.18) 



^2 ' - ^^2' ^^2' 

for s = ±1. The magnetization per unit cell is given by 



M = 



pi{k)dk 



p^i{k)dk, 



(2.19) 



where ps is the solution of Eq. (2.13) for each value of s. 

The equations ( ^.12| ) and (2.13) are solved numeri- 
cally. Figures 2(a-f) shows the magnetization curves for 
/3 = — 0.5, — 1 and —3. The open circles and open squares 
represent the results of high field and low field approxi- 
mations, respectively. 

§3. Quantum Monte Carlo Simulation 

Wc have also performed the quantum Monte Carlo 



simulation for the Hamiltonian (1.1) to check the re- 
gion of applicability of our boson approximation. The 
algorithm is the standard world line method using the 
checker board decomposition with local and global flip 
in the spatial and Trotter direction. The extrapolation 
with respect to the Trotter number is performed in the 
following way. For T = 0.1 and /3 = —3, for which the 
worst convergence is expected, we have taken the Trotter 
numbers L = 10, 12 and 18 and checked that the correc- 
tion term is almost proportional to 1/L^ as expected. 
The actual extrapolation is made by extrapolating the 
data for L = 12 and 18 linearly to The system 

size N is fixed to = 16 (32 spins). The measurement 
is done over typically 4 x 10^ Monte Carlo steps. The re- 
sults are shown in Figs. 2 (a-f) by filled circles. Typical 
error bars are smaller than the size of the symbols. 

§4. Results and Discussion 

The comparison of the data shows that the boson ap- 
proximation is fairly good for low field region < M<_0.1 
and near the saturation field 0.9<M < 1. In the inter- 
mediate regime, however, there exist considerable dis- 
crepancy between our approximation and numerical re- 
sults. One of the origin of the discrepancy is attributed 
to the finite band width of magnons. If we approximate 
them by the unbounded parabolic bands as Eqs. (2.11) 
and (2.16), we overestimate the number of magnons with 
large | g | in the limit of high (low) magnetic field within 
the low (high) field approximation. In the low field 
approximation, however, for the intermediate magnetic 
field, the number of magnons with g ~ tt are underesti- 
mated because the true dispersion is convex around this 



pointEP and these excitations are more easily excited than 
the parabolic case. For large negative /3, the convex part 
of the dispersion is absorbed in the two magnon contin- 
uumEP and the latter effect becomes less pronounced. It 
should be also noted that the multimagnon effect which 
is not correctly treated within the present approximation 
would also cause the discrepancy. For the high field ap- 
proximation, the — branch comes down as | /3 | decreases 
and it touches the + branch for f3 = —1. Therefore the 
contribution from the — branch is no more negligible 
for P ^ —1. This causes the underestimation of the 
magnon number for small /3 in the present single mode 
approximation. These arguments apply not only for the 
magnons in the ground state but also to the thermally 
excited magnons. Actually, the qualitative tendency of 
deviation does not depend on temperature as shown in 
Fig. 2. 

As an alternative method to map the spin chain onto 
the (5-function boson^as in the high field regime, Akutsu, 
Okunishi and HieidaEj^ have proposed to use the 2-body 
iS-matrix. They applied this method to the S — 1 
bilinear-biquadratic model and determined the interac- 
tion parameter c of the 5-function bose gas so that its 
2-body S'-matrix coincides with that of the spin chain 
in the long wave length limit. Surprisingly, their expres- 
sion for c turned out to be different from the one o\ 
tained by means of the Dyson-Maleev transformatior 
and diverges as the biquadratic term vanishes. Although 
the origin of the discrepancy is unclear, it may be re- 
lated with the order in which the long wavelength limit 
is taken and mapping onto the boson gas is made. Actu- 
ally, in the present model, the 2-body 5-matrix cannot 
be represented by a local phase shift and it is difficult 
to construct the 2-body S'-matrix before taking the long 
wave length limit. This problem is left for the future 
study. 

From the experimental side, several materials which 
realizes the spin- 1/2 ferromagnetic-antiferromagnetic 
alternating chain are synthesized. The magneti- 
zation curves are^ measured for [Cu(TIM)]CuCl4J$ 
(4-BzpipdH)CuCl3B' and methyl nitronyl nitroxide.LF 
The direct comparison of the present calculation and 
the experimental data is difficult, however, because of 
the presence of the interchain interaction which not only 
leads to the antiferromagnetic ordering at fairly high 
temperatures but also has influence on the magnetization 
curves in the paramagnetic phase. This kind of interac- 
tion cannot be treated within the present scheme which 
makes use of the Bethe ansatz solution. In this context, 
the synthesis of similar compounds with reduced inter- 
chain interaction is hoped for the direct comparison of 
the present calculation and experiment. On the other 
hand, the inclusion of other kinds of intrachain interac- 
tion is straightforward in the present scheme. Therefore 
our method gives a convenient way to calculate the low 
temperature magnetization curves of one-dimensional 
antiferromagnets in the low held and high held regime 
semi-quantitatively without resorting to elaborate nu- 
merical simulations. 

The author would like to thank S. Miyashita for pro- 
viding the source code of the quantum Monte Carlo sim- 
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ulation program. He also thanks M. Hagiwara for show- 
ing him his experimental data prior to publication and 
for stimulating discussion. Thanks are also due to Y. 
Natsume and T. Suzuki for sending their Monte Carlo 
data prior to publication and for invaluable comments. 
The numerical simulation is performed using the FA- 
COM VPP500 at the Supercomputer Center, Institute 
for Solid State Physics, University of Tokyo, the HITAC 
S820/80 and SR2201 at the Information Processing Cen- 
ter of Saitama University. This work is supported by the 
Grant-in- Aid for Scientific Research from the Ministry of 
Education, Science, Sports and Culture. 
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Fig. 1. The /3-dependence of A (open circles) and ^ (filled circles) 
calculated from the exact diagonalization data in ref ^. The solid 
lines are guides for the eye. 



Fig. 2. The magnetization curves of the present model. The re- 
sults of the Monte Carlo simulation, high field approximation 
and low field approximation are shown by the filled circles, open 
circles and open squares, respectively. The values of temperature 
T and fi are shown in the figure. 
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